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The probability distributions of discrete-time quantum walks have been often investi- 
gated, and many interesting properties of them have been discovered. The probability 
that the walker can be find at a position is defined by diagonal elements of the density 
matrix. On the other hand, although non-diagonal elements of the density matrices have 
an important role to quantify quantumness, they have not received attention in quantum 
walks. We focus on the non-diagonal elements of the density matrices for discrete-time 
quantum walks on the line and derive limit theorems for them. 
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1 Introduction 

The discrete-time quantum walk (QW) is considered as a quantum analogous of the random 
walk which expresses the random motion of particles, and many interesting properties of the 
QWs have been discovered [TJ HJ IH1 H] ■ Although QWs have simple dynamics in analogy with 
random walks, the behavior of quantum walkers is different from that of random walkers. 
The properties of the QWs are available to quantum information theory in which we can 
expect QWs to give a sizable contribution to quantum computer. Indeed, compared with 
classical search algorithms, the corresponding quantum algorithms based on QWs produce 
exponential speed-up. The difference between QWs and random walks is seen in probability 
distributions. The probability that the quantum walker on the line is found at position 
x € Z = {0, ±1, ±2, . . .} is defined by sum of diagonal elements in the density matrix at 
position x. The probability distributions of the QWs have been analyzed mainly and many 
limit theorems for them have been proved. Meanwhile, non-diagonal elements of the density 
matrices have got little attention in quantum walks and the behavior of them have not clarified. 
So, the aim of our paper is to get the characteristics of the non-diagonal elements in QWs. 
Non-diagonal elements of the density matrix, which are often called interference terms, show 
the distinction between quantum and classical states. If the particle takes mixed state, non- 
diagonal elements of the density matrix are necessarily null. That is, the interference terms of 
the density matrix for classical state are vanished. When the state of particles is not classical, 
they can become non-zero values. The non-diagonal elements of the density matrix have an 
important role when we characterize quantumness of particles. In quantum mechanics they 
are the essential values to express decoherence, which has been also investigated in QWs, 
both numerically and analytically [3j EJ |8] . 

The direct study for non-diagonal elements of the density matrices in QWs has not been 
done. Instead, we have analyzed the von Neumann entropy to quantify the entanglement, 
which is defined by density matrices [HI QUI EJ HH EH- Carneiro et al. [TT] numerically 
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computed the von Neumann entropy of QWs. In Abal et al. [9] and Annabestani et al. [10] . 
analytical results were obtained. Liu and Petulante [13] and Ide et al. [12] derived limit 
theorems for the von Neumann entropy. 

The present paper is organized as follows. In Sect. 2, we introduce the notations of 2-state 
QWs on the line and give the time evolution rule and probability distribution of the QWs. In 
Sect. 3, we present two limit theorems for interference terms of density matrices as our main 
result. Section 4 is devoted to proofs of the limit theorems. Summary is mentioned in the 
last section. In appendix A, we show limit theorems for 3-state Grover walks on the line. 

2 Definition of the 2-state QW on the line 

In this section we define the 2-state QW on the line. Let \x) (x G Z = {0,±1,±2, . . .}) be 
infinite components vectors which denote the position of the walker. Here, x-th component 
of \x) is 1 and the other is 0. Let \ipt(x)} G C 2 be the amplitude of the walker at position x at 
time t G {0, 1,2,.. .}, where C is the set of complex numbers. The walk at time t is expressed 
by 

|* t )=£|*)®hMar)). (1) 
The time evolution of the walk is described by the following unitary matrix: 



U 



u w li U 



(2) 



where Uj 1 j 2 G C(ji,j2 G {0,1}). In the present paper, we take the components of U as 
itoo = —uu = cos6>,uoi = Uio = sin 9 or i*oo = u n — cos 6*, — uoi — ifio = sin 6 with 
6 G (0, 2ir) (8 7r/2, tt, 3tt/2). Then the evolution is determined by 



\1> H .i{x)) = p\Mx + i)) + Q\Mz-i)) 

where P = |0) (0| U,Q = \1) (1| U and 



10) 



ID 



(3) 



(4) 



The probability that the quantum walker X t is at position x at time t, P(A t = x), is defined 
by 

i 

v{Xt=x)=Y,a\ptV)\j), (5) 

where pt{x) = \iftt(x)} (ipt(x)\ are the density matrices at position x at time t. 
3 Main results 

We will show two limit theorems for the non-diagonal elements of the density matrices pt{x). 
A limit theorem for the diagonal elements (j \ pt(x) \ j) (j G {0, 1}), that is probability distri- 
bution, has been given already by Konno [UJ [13] and Grimmett et al. [H] . Konno [UJ [TS] 
obtained the limit theorem for the probability distribution of X t /t as t — > oo by using path 
counting method. Grimmett et al. [16] also computed it by the Fourier analysis. From 
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their results, for the 2-state QWs starting from the origin with |-0o(O)) = a |0) + /3|1) and 
\a\ 2 + \f3\ 2 = 1, we see the following limit theorem. 
For uqoUoiUiqUu ^ and r = 0, 1, 2, . . ., we have 



lim^ff) £<j|M*) \j)= lim 



3=0 



Xt 

t 



y r g(y) dy, 



(0) 



where 



VI - Ko| 2 



tt(1 - x 2 )^\u 00 \ 2 - x' 



^(-luooMuool)^) 



X 1 



Ko| 2 



(7) 



Ia(x) = 1 ii x £ A, Ia{x) = x ^ A and E(Y) means the expected value of Y. We compare 
the limit function g{x) with probability distribution of the 2-state QW at time t = 1000 in 

Fig.ru 



Fig. 1. The limit function g(x) (thick line) and probability distribution of Xt/t at time t = 1000 
(thin line) with uqq = «oi = U 10 = — "11 = l/v^," = l/v2,/3 = — i/v2 



In the present paper, we focus on the non-diagonal elements (0| pt (a;) |1) in the density 
matrices pt(x), which are often called the interference terms in quantum mechanics, and derive 
limit theorems for them. At first, when we consider the 2-state quantum walks starting with 
any initial state, we can obtain the following relation between the interference terms and the 
moments of X t /t. 

Theorem 1 For r = 0, 1, 2, . . we have 

, lim E (jYm0\Pt(x)\l)) = (detCOf hm E 

t— >oo * — ' \ t / 2c t— >oo 

£CgZ 

where c = cos 9, s = sin 9 and $i(z) denotes the real part o/z£C. 

Moreover, if we assume that the 2-state quantum walks starting from the origin with |i/>o(0)) = 
a|0) + /3 |1) and \a\ 2 + \(3\ 2 = 1, the limit theorem is given as follows. 

Theorem 2 For r = 0, 1, 2, . . we see 

t ^E(f) r <°IM*)|i> = r y r {f {R) (y) + if {I) (y)} dy, (9) 
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where 



f (R) {x) =(det U) —x 



2c tt(1 - x 2 ) % /c 2_ 



{|a| 2 - |/3| 2 - (det (a/3 + a/3)} x] I(-\ c \,\c\)(x), 



{I) H3(a-/3) Vc 2 - x 2 
/ l W fZ^2 f (-|c|,|c|)(g), 



(10) 



(11) 



and 9(z) is the imaginary part of z 6 C. 

Figure [2] shows a comparison of the limit functions f^(x), f^(x) and non-diagonal elements 
of the density matrices pt(x). 



(a) fW(x),M({0\Pt(x)\l)) 



(b) /^J(x),^«0|p t (x)|l» 



Fig. 2. The limit function f (K) (x), f^{x) (thick line) and Sft((0] pt(«) S((0| pt(«) |1» at time 
t = 1000 (thin line) with uqo = -till = c = l/\/2,Uoi = u 10 = s = l/%/2,a = 1/^2, P = -i/V2 



4 Proof of Theorems 
4-1 Proof of Theorem^ 

Our proof is based on the Fourier analysis introduced by Grimmett et al. |16j . The Fourier 
transform \^ t (k)) (k € [— 7r,7r)) of \ipt(x)) is given by 



|# t (fc)) = 5>- fe hM*)>. 

By the inverse Fourier transform, we have 

^ dh 



(12) 



(13) 



From Eqs. ([3]) and (Tl"2"j) , the time evolution of |$t(fc)) becomes 

|*t+i(*0) = I7(k) \Mk)) , (14) 
where C7(fc) = R(k)U and fl(fc) = e lk |0) (0| + e~ lfc |1) (1|. Equation (H]) implies 

|*t(*)) = ^(*) t |*o(*)>. (15) 
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The eigenvalues Xj(k) (j £ {1, 2}) of the unitary matrix U(k) are 



icsin k — (— l) J 'yl — c 2 sin 2 fc (uqo = -«n = c, uqi = itio = s)> 



ccos fc — i(— l) J vi — c 2 cos 2 fc (uoo = Wn = c, — uoi = uio = s), 
and the eigenvectors corresponding to Xj(k) take the following form 



(16) 



-;ik 



-ccos 



fc-(-l)Vl-c 2 sin 2 fc 



se 



ik 



("00 = -"11 = C,Uqi = "10 = s), 



("00 = "11 = c, -"01 = "10 = s). 



(17) 



j {csinfc + (— l) J vl — c 2 cos 2 fc} 
By the eigenvalues \j(k) and the normalized eigenvectors \Vj{k)}, Eq. (|15[) can be rewritten 



as 



2 2 

\Mk)) = <«i(*)l*o(fc)) =X! A j( fc )' M fc )l*oW) . (18) 

From now, we will prove Theorem [T] At first, in order to express the left-side hand of Eq. ((SJ 
by the Fourier transform, we compute as follows: 

J2 ^((o| PtW |i» = i£* r {(o| |i) + <i| pt (x) |o)} 



x6 

1 

2 



(if t (k)\j+(D r \^ t (k)) 



dk 
2^ 



% r ( e wi («j(*)i*o(fe)) i 2 Mk)\j+\vj(k)) 



3 = 1 



(X 1 (k)X 2 (k)) t h 2 (k) r (»a(fc)|*o(A)) ("i(fc)| J + |" 2 (fc)) ("i(fc)|*o(fc)> 



+ (Ai(fc)A 2 (fc)) t /ii(fe)'" («i(*)|*o(*)) (" 2 (fc)| J+|"i(fc)) (" 2 (fc)|*o(fc)> 

+ o(t r " 1 ), 



dfc 
2^ 



(19) 



where J+ = |0) (1| + |1) (0|, h 3 (k) = DX 3 {k) / X^k) {j e {1,2}), L> = z(d/dfc) and (t) r = 
£(i — 1) x ■■• x (t — r + 1). By using the Riemann-Lebesgue lemma, we have 

( lim^g) r K((0|p t (x)|l»=i f J ^^(fcr|^(fc)|* (fc))| 2 (",(fc)|J + K(fc)) 1 



(20) 



Since we easily find (vj(k)\ J + \vj(k)) = (det U)-hj(k), 



f l_im^(|) r K(<0|p t (a:)|l>)=(det t/) 



=(det £/")— lim ] 

2c t-K» 



x;^(*r +i K«i(*)i*o(fc))i s 

r + ll 



2^ 



(21) 
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□ 

4.2 Proof of Theorem\M 

The real part f R \x) can be obtained from Theorem Q] and the result in Konno [TH [15] as 
follows: 



Hm^(|)>(0|Mz)|l>)^dett7), 



y 



r+l W 



1 



{|a| 2 - |/3| 2 - (detuf- (a/3 + a/3)} v] h-\c\M){y)dy 



y r (detU)-y 



2c tt(1 - y 2 )^Jc 2 - y 



2 



1 



{H 2 - |/3| 2 - (det[/)^(a/3 + a/?)}y] / ( _| c |,| c |)(»)di/. 

(22) 

Next, we calculate the imaginary part f^\x). In a similar fashion as the proof of Theorem[TJ 



°°xez J - ,r I 3=1 I 

(23) 



where denotes the imaginary part of z e C and J~ = |1) (0| — |0) (1|. Note that the 
initial state becomes |*o(*)) = l^o(0)) = a|0) + /?|1). After putting hj(k) =y(j € {1,2}), 
we get 



(f)^((0|^(-)|l>) = f I/ r]f! |^ • ^?A-W.H)(»)^ ( 24 ) 
Combining Eqs. (j2"2")l and we complete the proof of Theorem [3J 



□ 

5 Summary 

In this section, we mention discussion and conclusion for our results. In the present paper, 
we investigated non-diagonal elements in the density matrices of discrete-time 2-state QWs 
on the line. Non-diagonal elements of the density matrix are called interference terms and 
they are important values to quantify decoherence of quantum particles. In Sect. [3] we gave 
two limit theorems for the interference terms of density matrices at position x £ Z in 2- 
state QWs. One is the limit theorem for the QW starting with any initial state, the other 
is the limit theorem when the walker starts from the origin. Theorem [T] showed the relation 
between the interference terms of the QW starting with any initial state and the r-th moments 
(r = 0, 1, 2, . . .) of Xtjt. From Theorem (3] if the walker starts from the origin, the limit 
function which denotes the imaginary part of the interference terms always becomes an even 
function regardless of the initial state at the origin. One of the interesting problems is to 
get the relations between decoherent process and QWs, that is, expression of transition from 
quantum state to classical state by QWs. Chisaki et al. [18] studied crossover from QW 
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to classical random walk. So, it would be interesting to discuss the relation between our 
results and their results. If we can construct the decoherent process by QWs, innovative 
investigations may start in quantum mechanics. 
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Appendix A 

In this appendix, we consider the 3-state Grover walk on the line and will show two limit 
theorems for the non-diagonal elements of the density matrices. The difference between the 
2-state QWs defined in Sect. [5] and the 3-state Grover walk is whether localization occurs 
or not. The probability distribution of the 2-state QWs doesn't localized, while the 3-state 
Grover walk is known as one of the models in which localization can occur. The 3-state Grover 
walk on the line is defined as follows. The amplitudes of the walker at position x G Z at time 
t G {0, 1,2,.. .} are given by \ipt(x)) G C 3 . The time evolution of the walk is determined by 



the Grover matrix G 



E2 ^"^^ 
ji =o j 2 



MiM-EU li) 01 with 





" 1 " 




" " 




" " 


10} = 





,11} = 


1 


,|2) = 

















1 



(A.l) 

1 

The amplitudes \^t(x)) evolve by the following rule: 

\^t+i(x)) = P \Mx + 1)} + Pi \Mx)) + Pi U't{x - 1)} , (A.2) 

where P = |0}(0|G,Pi = |1}(1|G,P 2 = |2)(2|G. Inui et al. HZ] focused on the prob- 
ability distribution P(X t — x) — Ylj=o 01 Pt{ x ) \j) an d got a long-time limit theorem for 
Exez(!) r Ei=oO'IMaO|j) = E[(*t/t) p ] (see Eq. (16) in [17). Figure EU shows proba- 
bility distribution of the 3-state Grover walk at time t — 150 with |-0o(O)) = l/v3 |0) + 
i/V3 \l) +i/V3 |2). 



-150 -75 75 150 



Fig. A.l. Probability distribution of the 3-state Grover walk at time t 
origin with the initial state |f/>o(0)} = l/VS |0> + i/V3 |1> + i/V3 |2) 



150 starting from the 



On the other hand, we derive limit theorems for the interference terms \ pt{x) | j2 ) ((ji, J2) G 
{(0, 1), (0, 2), (1, 2)}). The proofs of theorems are given as well as the ones of Theorems Q] and 
[21 So, we will omit them. 

For the 3-state Grover walk starting with any initial state, we can obtain the relations 
between the interference terms and the moments of Xt/t. 

Theorem A.l As t — > 00, the summations of (x/t) r $l((ji \ pt(x) | j'2) ) (r — 0,1,2,...) all 
over ifZ are calculated as follows: 

1. For r = 0, we have 



lim J2mn\Pt(x) |j 2 }) 



a(R) 
^01 



lim^oo E - \ lim t _ 



E 



A (R) , 1 5 1- 



lim* 



E 



lim* 



(ii = o,j 2 = i), 

(.71 = 0,72=2), 
(ji = l,i 2 =2), 



(A3) 
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where 



3132 



{ r-^\\Hk)\Mk))\ 2 {v{k)\J+ j2 \v{k)) (0'1,J2) = (0,1),(1,2)), 

f {*IM*OI*o(*0>l>(*)l^>(*0> 

-i|( w (fc)|* (fc))| 2 } (0'i,J2) = (0,2)), 



(A.4) 

and |* (fc)) = E^e - '** IM^)> and J+ J2 - (j 2 | + |j 2 ) 77ie sector \v{k)) is 
the normalized eigenvector corresponding to the eigenvalue A = 1 of the unitary matrix 
R(k)G with R(k) = e ik |0) (0| + |1) (1| + e~ ik |2) (2|: 



\v(k)) 

2. For r = 1, 2, . . we have 

t !™E(f) rR (0'i|p*(») |J2» 



5 + cos fc 



(l + e- jfe )/2 



= < 



(A.5) 



| lim^oo E 


> ) r+r - 




E [( 




(ii 


= 


J2 


= 1) 


i linit^co E 


l(^) r " 


5 
4 


lim^oo E 




y+2~ 




(ii 


= 


J2 


= 2) 


-ilim t ^ 00 E[(^) 


r+1" 


- \ linit^oo E 


\¥) 


r+2" 


(ii 


= i 


J2 


= 2) 



(A.6) 



(A.7) 



Particularly, when the 3-state Grover walk starts from the origin with IV'o(O)) = a |0) + /3 |1) + 
7 |2) and \a\ 2 + \fl\ 2 + \j\ 2 = 1, we can get the specific limit theorem for the interference terms. 
Theorem A. 2 For r = 0, 1, 2, . . ., we see 

^E(f) r 0ilft(*)|j»> = /~ y r {fS 2 (y) + ^!l(y)} d ^ ( A - 8 ) 



where 



/;» = 



l A 12 ^o(*) - ^Vr^ I^-*,*)^ 0i = U 2 = 2), 
A$8 (x) (Wx)vCT ^ J( __^^ )(a;) (ji = = 2)) 



(A.10) 
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and Sq{x) denotes Dirac's S -function at the origin. The coefficients are computed as follows: 



and 



A 



(i) 



s/6 
Mi 



7+1 



+ (l - 2 ^f) 3? {(2a + /3)(27 + /3)} (ji = 0,j a = 1), 
-#(h+ff + |7+f| 2 ) 

+ (2 - ii^l) n{(2a + (3) (27 + £)} tfi = 0, j 2 = 2), 



36 



7+f 



+ (l - 2^1) 3? {(2a + /3)(2 7 + 0)} (ji = hh = 2), 



2 - s^eA 3 



J 3(a/3 + /? 7 + 2 7 a) (ii = 0,j 2 = 1), 



4 _ 
2 



3(a/? + /3 7 + 2 7 a) (j, = 0, j 2 - 2), 
)3(a/? + ^ 7 + 2 7 a) (ii = l,j 2 = 2), 



5\/S\ 
6 



c = |a + 7 | 2 + 2|/3| 2 , 

Cl = 2(-|a~/3| 2 + | 7 -/3| 2 ), 

c 2 = |a- 7 | 2 -23?{(2a + /3)(2 7 +^))}, 

d -3{(a + 7 )^}, 

di - 3{a(/3 + 7 ) + (a + /3) 7 } . 



(A.H) 



(A.12) 



(A.13) 
(A.14) 
(A.15) 
(A.16) 
(A.17) 



We should note = A 12 , = A^ 2 ■ When Theorems I A. II and IA.2I are proved, we con- 
sider J+j = \n) (j2\ + \h) (ill (resp. Jr.^ = \j 2 ) O2I) instead of J+ = |0) (1| + |1) (0| 

(resp. J~ = |1) (0| — |0) (1|) in the proofs of Theorems [TJ and [H Examples of Theorem IA. 21 
are presented in Figs. IA.2I and I A. 31 
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(a) (ji,j 2 ) = (0,1) (b) Cn,j 2 ) = (0,2) (c) (j u j 2 ) = (1,2) 



Fig. A.2. The limit function f\ ! (x) (thick line) and 5R((ji | pt(a;) |j 2 >) at time t = 1000 (thin line) 
with a = 1/V3, P = 7 = i/>/3 




(a) 0'i.Ja) = (0,1) (b) (ji, j 2 ) = (0,2) (c) (ji,j 2 ) = (1,2) 



Fig. A.3. The limit function (a) (thick line) and S>(0'l I Pt 0*0 |j2» at time t = 1000 (thin line) 
with a = 1/V3, 13 = -/ = i/V3 



